Electron spin relaxation in paramagnetic Ga(Mn)As quantum wells 
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Electron spin relaxation in paramagnetic Ga(Mn)As quantum wells is studied via the fully micro- 
scopic kinetic spin Bloch equation approach where all the scatterings, such as the electron-impurity, 
electron-phonon, electron-electron Coulomb, electron-hole Coulomb, electron-hole exchange (the 
Bir-Aronov-Pikus mechanism) and the s-d exchange scatterings, are explicitly included. The Elliot- 
Yafet mechanism is also incorporated. From this approach, we study the spin relaxation in both 
n-type and p-type Ga(Mn)As quantum wells. For n-type Ga(Mn)As quantum wells where most 
Mn ions take the interstitial positions, we find that the spin relaxation is always dominated by 
the DP mechanism in the metallic region. Interestingly, the Mn concentration dependence of the 
spin relaxation time is nonmonotonic and exhibits a peak. This is due to the fact that the mo- 
mentum scattering and the inhomogeneous broadening have different density dependences in the 
non-degenerate and degenerate regimes. For p-type Ga(Mn)As quantum wells, we find that the Mn 
concentration dependence of the spin relaxation time is also nonmonotonic and shows a peak. Dif- 
ferently, the cause of this behaviour is that the s-d exchange scattering (or the Bir-Aronov-Pikus) 
mechanism dominates the spin relaxation in the high Mn concentration regime at low (or high) 
temperature, whereas the DP mechanism determines the spin relaxation in the low Mn concentra- 
tion regime. The Elliot- Yafet mechanism also contributes to the spin relaxation at intermediate 
temperatures. The spin relaxation time due to the DP mechanism increases with increasing Mn 
concentration due to motional narrowing, whereas those due to the spin-flip mechanisms decrease 
with it, which thus leads to the formation of the peak. The temperature, photo-excitation density 
and magnetic field dependences of the spin relaxation time in p-type Ga(Mn)As quantum wells are 
investigated systematically with the underlying physics revealed. Our results are consistent with 
the recent experimental findings. 

PACS numbers: 72.25.Rb, 75.50.Pp, 71.10.-w, 71.70.Ej 



I. INTRODUCTION 

Semiconductors doped with magnetic impurities have 
intrigued much interest since the invention of ferro- 
magnetic III-V semiconductors due to the possibility 
of integrating both the magnetic (spin) and charge 
degree of freedom on one chip~£&££&2- Many new 
device conceptions and functionalities based on these 
materials are proposed, and the material properties 
together with the underlying physics are extensively 
Studied4.^ii0 i 1 i 1 i 2 . i 1 3ii4 i i5 1 i6 Specifically, ferromagnetic 

Gai_ a; Mn 2 ;As has been used as a highly efficient source 
to inject spin polarization into GaAsii and magnetic tun- 
neling junctions based on ferromagnetic Gai-^Mn^As 
can achieve very high magnetoresistanceJ^ Besides, 
the ability to detect the magnetic moment via Hall 
measurement o 4 ' 19 and to control it via gate- voltage^ and 
laser radiation 2 ^ opens the way for incorporating opto- 
electronics with magnetism. Magneto-optical measure- 
ments, which could characterize the spin splitting of car- 
riers due to both the external magnetic field and the s- 
d or p-d exchange field, provide important information 
about the microscopic properties of the carriers, such as, 
the <?-factor, the s-d and p-d exchange coupling constants, 
as well as the electron spin relaxation time (SRT). Such 
measurements have recently been performed in dilutely- 



doped paramagnetic Ga(Mn)As quantum wellsi 22 i 23 ' 24 i 25 
Although many properties of Ga(Mn)As have been ex- 
tensively studied, the electron spin relaxation has not 
yet been well understood even in dilutely doped param- 
agnetic phase. This is the aim of this investigation. We 
focus on (001) Ga(Mn)As quantum wells. 

Electron spin relaxation in non-magnetic GaAs has 
been extensively studied and three main spin relaxation 
mechanisms have been recorgnized for decades^ the 
D'yakonov-Perel' (DP) mechanism^ the Bir-Aronov- 
Pikus (BAP) mechanism 2 ^ and the Elliot- Yafet (EY) 
mechanism. 29 Usually, the DP mechanism dominates the 
spin relaxation in n-type quantum wellsi 10 ' 30 The BAP 
mechanism was believed to be most important at low 
temperature in intrinsic and p-type quantum wells for a 
long time. 10 ' 26 Recently, Zhou and Wu showed that the 
BAP mechanism was exaggerated in the low temperature 
regime in previous treatments based on the elastic scat- 
tering approximation, where the Pauli blocking was not 
considered^ It was then found that the BAP mechanism 
is less efficient than the DP mechanism in intrinsic quan- 
tum wells and p-type quantum wells with high photo- 
excitation density. A similar conclusion was also ob- 
tained in bulk GaAs very recently.— Previously, the EY 
mechanism was believed to dominate the spin relaxation 
in heavily doped samples at low temperature. However, it 
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was shown to be unimportant in bulk GaAs by our recent 
investigation^ Whether this is still true in quantum-well 
systems remains unchecked. Moreover, in paramagnetic 
Ga(Mn)As quantum wells, things are more complicated: 
(i) All the three mechanisms could be important as the 
material is heavily-doped with Mn and the hole density is 
generally very highj^ 2 - (ii) An additional spin relaxation 
mechanism due to the exchange coupling of the electrons 
and the localized Mn spins (the s-d exchange scattering 
mechanism) may also be important £ In this work, we will 
compare different spin relaxation mechanisms for various 
Mn concentrations, temperatures, photo-excitation den- 
sities and magnetic fields. 

In Ga(Mn)As, the Mn dopants can be either substi- 
tutional or interstitial: the substitutional Mn accepts 
one electron, whereas the interstitial Mn releases two. 
Direct doping in low-temperature molecular-beam epi- 
taxy growth gives rise to more substitutional Mn ions 
than interstitial ones, which makes the Ga(Mn)As a p- 
type semiconductor ! 8 ' 9 ' 22 Recently, it was found that in 
GaAs quantum wells near a Ga(Mn)As layer, the Mn 
dopants can diffuse into the GaAs quantum well, where 
the Mn ions mainly take the interstitial positions, mak- 
ing the quantum well n-type i 23 ' 24 ' 33 The experimental 
results also show interesting features of the SRT. 

We apply the fully microscopic kinetic spin Bloch 
equation (KSBE) approac h 34 ' 35 ' 36 to investigate the 
spin relaxation in paramagnetic Ga(Mn)As quantum 
wells. The KSBE approach explicitly includes all 
relevant scatterings, such as, the electron-impurity, 
electron-phonon, electron-electron Coulomb, electron- 
hole Coulomb, electron-hole exchange (the BAP 
mechanism) and s-d exchange scatterings. Previously, 
the KSBE approach has been applied to study the 
spin dynamics in semiconductor and its nanostructures 
where good agreements with experiments have been 
achieved and many predictions have been confirmed by 

c:: p Cr - mcn ^ :: 30.31.32.34.35.36.37.38.39.40.41.42.43.44.45.46.47.48.49 

In this work, we apply the KSBE approach to both 
n- and p-type paramagnetic Ga(Mn)As quantum wells 
to study the electron spin relaxation. We distinguish 
the dominant spin relaxation mechanisms in different 
regimes and our results are consistent with the recent 
experimental findings. 

This paper is organized as follows: In Sec. II, we set 
up the model and establish the KSBEs. In Sec. Ill 
we present our results and discussions. We conclude in 
Sec. IV. 



II. MODEL AND KSBES 

We start our investigation from a paramagnetic [001] 
grown Ga(Mn) As quantum well of width a in the growth 
direction (the z-axis). A moderate magnetic field B is 
applied along the :r-axis (the Voigt configuration). It is 
assumed that the well width is small enough so that only 
the lowest subband of electron and the lowest two sub- 



bands of heavy-hole are relevant for the electron and hole 
densities in our investigation. The barrier layer is chosen 
to be Alo.4Gao.6As where the barrier heights of electron 
and hole are 328 and 177 meV respectively^ The enve- 
lope functions of the relevant subbands are calculated via 
the finite-well-depth mode h 30 ' 31 

The KSBEs can be constructed via the nonequilibrium 
Green function method 51 and read 



dtpk = <9t/5k| 



<9t/S k | 



(1) 



with /5k representing the single-particle density matrix 
whose diagonal and off-diagonal elements describe the 
electron distribution functions and the spin coherence 
respectively^ 4 - The coherent term is given by (fi = 1 
throughout this paper) 

dtPklcoh = 



(g e ^ B B 



h(k)).| 



(k),Pk 



(2) 



in which [ , ] is the commutator. g e is the electron g- 
factor. h(k) represents the spin-orbit coupling (SOC), 
which is composed of the Dresselhaus 5 - 2 - and Rashba 53 
terms. In symmetric GaAs quantum well with small well 
width, the Dresselhaus term is dominant 5 - 4 - and 

h(k) = 7D (k x [ky - <fc 2 )) , k y «A£) - k 2 x ) , 0) . (3) 

Here {k 2 ) represents the average of the operator 
— (d/dz) 2 over the state of the lowest electron sub- 
band and 7d = 11-4 eV-A 3 is the Dresselhaus SOC 
coefficient^ The mean-field contribution of the s-d ex- 
change interaction is given by 



H% = -N Mn a(S) ■ j ' 



(4) 



where (S) is the average spin polarization of Mn ions and 
a is the s-d exchange coupling constant. For simplicity, 
we assume that the Mn ions are uniformly distributed 
within and around the Ga(Mn)As quantum well with a 
bulk density N Mn . Shf (k) = ~IZq, 9 , v i,q m \I{nz)? Pk- q 
is the Coulomb Hartrce-Fock (HF) term, where I(iq z ) = 
J dz\t; e (z)\ 2 e tq * z is the form factor with £ e (z) standing for 
the envelope function of the lowest electron subband. 35 
Vq t q z is the screened Coulomb potential. In this work, we 
take into account the screening from both electrons and 
holes within the random phase approximation! 3 - 1 - 

The scattering term <9t/3k| scat consists of the electron- 
impurity, electron-electron Coulomb, electron-phonon, 
electron-hole Coulomb, electron-hole exchange and s- 
d exchange scatterings. The expressions of all these 
terms except the s-d exchange scattering can be found in 
Ref.l3ll However, the expression of the electron- impurity 
scattering term with the EY mechanism included has not 
been given in that paper, which we will present later in 
this paper. The s-d exchange scattering term is given by 
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-irN Mn a 2 I s C?Mn (-??! ~m) 5 ( £ k - £ k') 

771772 k' 

s m p^,s m p^ - s^p<s m pl + H.c. 



(5) 
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Here p> = I - p k , 
±Tr(5^p Mn ) and I„ 



ft 



S 11 and 



(7? = 0, ±1) are the spin ladder operators with 5° = S z , 
5 ± = S x ± iS y , s° — 2s z and 



. ;± = s x ± is y . pMn is the 
Mn spin density matrix. £ k = k 2 /2m* is the electron 
kinetic energy with to* denoting the effective mass. The 
equation-of-motion for Mn spin density matrix consists 
of three parts 5 t p Mn = d t pMn\ coh + d t pMn\ scat + dtf>Mn\ rcl - 
The first part describes the coherent preces- 
sion around the external mangetic field and 
the s-d exchange mean field, $tPMn| co h = 

-* ffMnMBB • S - aJ2k Tr (f-Pk) • S, PMn • The 
second part represents the s-d exchange scattering 
with electrons d t p M „| scat = £„,„, k <K £ k - 



Ev)Tr(s-^p<s-vip>) (S^S^pun 



77iT/ 2 k 
S^PMnS^) 



H.c. . The third term characterizes the Mn spin 

relaxation due to other mechanisms, such as, p-d 
exchange interaction with holes or Mn-spin-latticc 
interactions, with a relaxation time approximation, 
<9t/°Mn| rc i = - {mn - pM n ) / T Mn - Here represents 
the equilibrium Mn spin density matrix, rjvin is the Mn 
spin relaxation time, which is typically 0.1^10 ns^ 
In our calculation we take TMn = 1 ns. At t = 0, the 
Mn spin density matrix is chosen to be the equilibrium 
one f>1i n . The Mn spins can be dynamically polarized 
via the s-d exchange interaction, and feedback to the 
electron spin dynamics. However, we find that this 
process affects the electron spin dynamics little. Hence, 
the choice of rjvin does not affect our discussions on 
electron spin dynamics. 

The s-d exchange scattering r st j can be obtained ana- 
lytically. In the absence of an external magnetic field, 
the spin polarization of the electron system is always 
along the z-direction. As the s-d exchange interaction 
conserves the spin polarization of the total system, the 
spin polarization of Mn ions (which is assumed to be zero 
initially) can only be along the z-direction. By keeping 
only the diagonal term of pMni from the Fermi Golden 
rule, the spin relaxation time due to the s-d exchange 
scattering can be obtained directly: 



r«,i = \ - N Mn I s a 2 m* [5(5+1) - (S 2 Z )} 
12 



35NMnCt 2 m*I s ' 



(6) 



where 5 = 5/2 is the spin of the Mn ion. It is evident that 
r s d is independent of temperature and electron density, 
but is inverse proportional to Mn concentration iVMn, 
the square of the s-d exchange coupling a 2 and I s — 
j dz\^ e (z)\ A . I s is determined by the confinement of the 
quantum well. For an infinite-depth- well I s = thus 
T s( j is proportional to the well width a. 

After incorporating the EY mechanism, besides the or- 
dinary spin-conserving term, there are spin-flip terms. 
For electron-impurity scattering these additional terms 



dtPk 



EY 



k' 



(1) A > A« 



n.c. , (7) 



where n, = N^ n + 4A^ n + n i0 with N^ n , A^ n , n i0 
being the densities of substitutional Mn, interstitial 
Mn and non-magnetic impurities, respectively due to 



different charges. u£2 k , = f EJk-^ 5 J^fe)l 2 ? 
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and U, 



(2) 



Here A c 
/-.',. and Ago are the 



3m c B B (l-r,/3) Wlth ^ _ A so +E 

band-gap and the spin-orbit splitting of the valence band, 
respectively.— The spin-flip matrices are given by A k V k = 
[(k + k', 0) x &] x and A k 2 k , = [(k, 0) x (k', 0)] • &. It is 



noted that A k V k and A^,'^ contribute to the out-of-plane 
and in-plane spin relaxations, respectively. They are gen- 
erally different and therefore the spin relaxation due to 
the EY mechanism in quantum wells is anisotropic. The 
EY mechanism can be incorporated into other scatter- 
ings similarly. 3 - However, we find that the EY mecha- 
nism can be important only when the impurity density is 
high, where the electron-impurity scattering is most im- 
portant. Therefore, for simplicity, we include only the EY 
spin-flip processes associated with the electron-impurity 
scattering. 
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III. RESULTS AND DISCUSSIONS 

By solving the KSBEs numerically, we obtain the tem- 
poral evolution of the single particle density matrix p k 
and then the spin polarization along the z-axis, i.e., s z . 
From the decay of s z , the SRT is extracted. The initial 
spin polarization is chosen to be P = 4 %. The well width 
a = 10 nm. The external magnetic field B is zero unless 
otherwise specified. We use x to denote the Mn density, 
where N Mn = xN with N = O -1 = 2.22 x 10 22 cm" 3 
(fi is the volume of the unit cell in GaAs). The other 
material parameters used are listed in Table [J 56 i 57 i 58 

The value of the s-d exchange coupling in III(Mn)V 
materials is still in dispute. In bulk Ga(Mn)As, first- 
principle calculation gives the value Noa ~ 0.25 eV££ 
However, the experimental measurements in Ref . I22I show 
that Noa is in the range of [—0.21, —0.07] eV varying with 
quantum well width. In this paper, we choose Noa — 
—0.15 eV unless otherwise specified. 



Spin relaxation in n-type Ga(Mn)As quantum 
wells 



In this subsection, we study the electron spin relax- 
ation in n-type Ga(Mn)As quantum wells, where the Mn 
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TABLE I: Material parameters used in the calculation 



N Mn (10 n cm- 2 ) 



5.31 x 10 3 kg/m 3 ei4 1.41 x 10 9 V/r 



v a i 5.29 x 10 3 m/s 
ojlo 35.4 meV 



KO 12.9 Koo 10.8 

D 

Vst 2.48 x 10 3 m/s 

3 8.5 eV 

A so 0.341 eV E g 1.55 eV 

AElt 0.08 meV a 146.1 A 

g e -0.44 m* 0.067m 

gMn 2 g 5/2 




dopants mainly take interstitial positions. For simplic- 
ity, we neglect the substitutional Mil's and assume that 
electrons from Mn donors are all free electrons. We will 
discuss the situations that the quantum wells are either 
undoped or n-doped before Mn-doping. 

For quantum wells which are undoped before Mn- 
doping, N e = Y e Mn + Y ox where Y e Mn i s the density of 
electrons from Mn donors and N ox is the photo-excitation 
density. We choose iV cx = 10 10 cm~ 2 which is usually 
smaller than iV e Mn . The SRTs due to various mecha- 
nisms are plotted as function of x in Fig. []Ja). N a is 
chosen to be —0.25 eV, which is smaller (i.e., the s-d ex- 
change interaction is stronger) than the value measured 
by experiments^ However, even for such a strong ex- 
change coupling, the spin relaxation due to the s-d ex- 
change scattering mechanism is still much weaker than 
that due to the DP mechanism. It is further seen from 
Fig. QJa) that the BAP and EY mechanisms are also 
unimportant. Therefore, the SRT is determined by the 
DP mechanism. Interestingly, the SRT due to the DP 
mechanism tdp first increases then decreases with in- 
creasing x. The t-x curve thus has a peak. The underly- 
ing physics is that the SRT has different x (density) de- 
pendence in the non-degenerate and degenerate regimes. 
Similar behavior has been found in bulk non-magnetic 
III- V semiconductors in Ref. [32j very recently. Let us first 
recall the widely used expression, tdp = l/[(|h(k)| 2 )T p ] 
((...} denotes the ensemble average), which is derived 
within the elastic scattering approximation and is only 
correct qualitatively^ The expression contains two key 
factors of the DP spin relaxation: (i) the inhomoge- 
neous broadening from the k-dependent spin-orbit field 
~ (|h(k)| 2 ); (ii) the momentum scattering time t p . The 
SRT due to the DP mechanism increases with increasing 
momentum scattering but decreases with increasing inho- 
mogeneous broadening. It should be mentioned that for 
this system, N e w n.i/2 = 2xNq (Note that the charge 
number Z of the Mn ion is included in rij as Z 2 . For 
interstitial Mn, which acts as a double donor, Z = 2.) 
In the small x (low density) regime, the electron system 
is in the non-degenerate regime, and the distribution is 
close to the Boltzmann distribution. Therefore, the in- 
homogeneous broadening of the k dependent spin-orbit 
field ~ (|h(k)| 2 ) changes little with electron density N e 
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FIG. 1: (Color online) SRT r due to various mechanisms in n- 
type Ga(Mn)As quantum wells which are (a) undoped or (b) 
n-doped before Mn-doping as function of Mn concentration x 
at 30 K (•) and 200 K (□). Red solid curves: the SRT due to 
the DP mechanism tdp; Green dotted curves: the SRT due 
to the EY mechanism tey; Brown dashed curves: the SRT 
due to the BAP mechanism tbap; Blue chain curve: the SRT 
due to the s-d exchange scattering mechanism r s d . The Fermi 
temperature of electrons T§ is ploted as black curve with A 
(the scale of Tp is on the right hand side of the frame) and 
T§ =T for both T = 30 and 200 K cases are plotted as black 
dashed curves. We also plot the scale of the electron density 
from Mn donors N^ 111 on the top of the frame. 



(hence x). On the other hand, in the non-degenerate 
regime the electron-electron scattering increases with in- 
creasing electron densit y 60 ' 61 (thus x). Moreover, the 
electron-impurity scattering also increases with increas- 
ing x as the impurity density increases. Therefore, tdp 
increases with increasing x (motional narrowing). In 
large x (high density) regime, the electron system is in 
the degenerate regime, the inhomogeneous broadening 
changes as kp oc N e oc x (fcf oc cx a: 3 ) if the lin- 
ear (cubic) term dominates the SOC. On the other hand, 
the electron-electron scattering decreases with increasing 
electron density (thus x) in the degenerate regimei 60 ' 61 
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Besides, the electron-impurity scattering increases slower 
than Ni oc x because the scattering cross section de- 
creases as the electron (Fermi) energy increases. Thus 
for both the linear- and cubic-term dominant cases, tdp 
decreases with increasing x in the large x regime. Conse- 
quently, tdp first increases then decreases with increas- 
ing x and a peak is formed in the crossover regime where 
T Tp (Tp is the electron Fermi temperature). It is 
seen from Fig. [Ha) that for both T = 30 and 200 K 
cases, the peaks indeed appear at T ~ Tp. It should be 
pointed out that the situation here is different from that 
in Ref . H2, where the density dependence of the SRT also 
has a peak in intrinsic quantum wells at room temper- 
ature. In that case, the impurity density is rather low 
and the relevant momentum scatterings are the carrier- 
carrier Coulomb and electron-phonon scatterings. In the 
situation here, the impurity density is extremely high 
(jii = 2N e ) and the relevant momentum scattering is the 
electron-impurity scattering. 

We now turn to the situation that the quantum wells 
are n-doped before Mn-doping. In this case, N e = 
Ni + A e Mn + A ox , where N l e denotes the density of 
the electrons from other dopants which is chosen to be 
10 11 cm~ 2 . We assume that the other dopants are far 
away from the quantum wells, so that they contribute 
little to the electron-impurity scattering, corresponding 
to the genuine case of modulation doping. However, the 
Mn ions are doped in the quantum wells i 23 ' 24 The photo- 
excitation density is N cx — 10 10 cm~ 2 . The results are 
plotted in Fig. QJb). As the density of the photo-excited 
holes is much smaller than the electron density, the BAP 
mechanism is obviously negligible and thus not plotted in 
the figure. From the figure, it is noted that the EY and s- 
d exchange scattering mechanisms are also insignificant. 
Consequently, the spin relaxation is still dominated by 
the DP mechanism. Similar to that in Fig. QJa), the 
SRT due to the DP mechanism tdp first increases then 
decreases with increasing x. For the case of T = 200 K, 
the peak of the SRT is still around T — Tp. However, for 
the case of T = 30 K, the peak moves to a larger x value 
compared to that in Fig. HJa). This can be understood 
by noting that the electrons have two sources: the Mn 
donors and other dopants. For the case of T = 200 K, the 
crossover of the non-degenerate and degenerate regimes 
takes place around T ~ Tp, where the corresponding x 
is 10~ 5 . At such x, electrons are mainly from the Mn 
ions rather than from other dopants. Thus the situation 
is the same as that in Fig. [IJa) and the peak appears at 
T ~ Tp. However, in the case of T = 30 K, for all x in the 
figure, T| is larger than T and the situation is hence dif- 
ferent. The t-x behavior in this case can be understood 
as follows: For x < 10~ 6 , electrons are mainly from the 
other dopants and N e changes slowly with x, thus the 
inhomogeneous broadening varies slowly with x. On the 
other hand, the electron-impurity scattering increases as 
the impurity density increases n, ~ Ax No. At low tem- 
perature (T < Tp), the electron-impurity scattering can 
be important even when rii < A^. 30,32 Therefore, the mo- 



mentum scattering increases with x significantly. Conse- 
quently, tdp increases with increasing x. For x > 10~ 5 , 
electrons mainly come from the Mn donors. The scenario 
becomes the same as that in the case of Fig. [Ha) and the 
SRT decreases with increasing x as Tp is much larger 
than T. Consequently, the peak is formed in the range 
10~ 6 < x < 10 -5 at x — 3 x 10~ 6 , which is larger than 
that in the case of T = 30 K in Fig. HJa). 

It should be mentioned that our results are consistent 
with the latest experimental finding that in the low Mn 
concentration regime the SRT increases with increasing 

xM& 



B. Electron spin relaxation in p-type Ga(Mn)As 
quantum wells 

In this subsection, we discuss the electron spin relax- 
ation in p-type Ga(Mn)As quantum wells. Both sub- 
stitutional and interstitial Mn ions exist in the system. 
Each substitutional Mn donates one hole, whereas each 
interstitial Mn compensates two holes. For simplicity, we 
assume that all the holes are free. The ratio of the hole 
density Nh to the Mn density Am„ is obtained by fit- 
ting the experimental data in Ref. [22], as shown in Fig. [2] 
From these densities, according to charge neutrality, the 
densities of substitutional Mn N^ n and interstitial Mn 
N^ An are determined. The photo-excitation density is 
chosen to be A cx = 5x 10 10 cm~ 2 unless otherwise spec- 
ified. 




0.02 0.04 0.06 0.08 0.1 0.12 
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FIG. 2: (Color online) Ratio of the hole density to the 
Mn density N^/Nmu vs. the Mn concentration x in p-type 
Ga(Mn)As quantum wells. The black dots represent the ex- 
perimental data. The red solid curve is the fitted one. The 
hole density Nh is also plotted (the blue dashed curve). Note 
that the scale of Nh is on the right hand side of the frame. 
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FIG. 3: (Color online) SRT r due to various mechanisms and 
the total SRT in p-type Ga(Mn)As against the Mn concen- 
tration x at (a) T = 5, (b) 50 and (c) 200 K. We also plot the 
scale of JVivin on the top of the frame. 



1. Mn concentration dependence of the SRT 

We first study the Mn concentration dependence of the 
SRT. In Fig. [3J we plot the SRTs due to various mech- 
anisms and the total SRT as function of the Mn con- 



centration x at T — 5, 50 and 200 K. It is noted that 
the total SRT first increases and then decreases with in- 
creasing x and there is a peak at x ~ 3 x 10 -5 . Remark- 
ably, the spin relaxation at large x is not dominated by 
the DP mechanism, but by the s-d exchange scattering 
(or the BAP) mechanism at low (or high) temperature. 
At medium temperature, the EY mechanism also con- 
tributes for large x. The SRT due to the DP mechanism 
increases with increasing x, whereas those due to the s-d 
exchange, EY and BAP mechanisms decrease with in- 
creasing x. Consequently a peak is formed. It should 
be pointed out that the underlying physics here is differ- 
ent from that in the case of n-type Ga(Mn)As quantum 
well where the peak is solely due to the electron den- 
sity dependence of the DP spin relaxation. It should be 
mentioned that the peak position is x ~ 10 -4 , which is 
consistent with that observed in Ref. I22I 

Let us now turn to the x dependence of the SRT due 
to various mechanisms. The increase of tdp with in- 
creasing x is due to the increase of the electron-impurity 
and electron- hole scatterings (motional narrowing). For 
the EY mechanism, the SRT decreases as the spin-flip 
scattering increases with increasing impurity density [see 
Eq. 0] . The s-d exchange scattering increases with in- 
creasing x as the Mn density increases [see Eq. (JS])]. The 
x dependence of the the BAP spin relaxation is more 
complicated. To facilitate the understanding, we plot 
tbap from the full calculation and that from the calcu- 
lations without the Pauli blocking of electrons (holes) 
in Fig. [H It is seen that for x < 10~ 6 , tbap changes 
little with x. This is due to the fact that the holes 
from Mn dopants are much fewer than those from photo- 
excitation, and hence Nh changes little with x. So does 
tbap- For larger x, tbap first decreases then increases a 
little and finally saturates with increasing x at T = 5 K. 
It is noted that without the Pauli blocking of holes, tbap 
decreases with increasing x rapidly, which indicates that 
the slowdown of the decrease and the saturation of tbap 
are due to the Pauli blocking of holes. It is further shown 
that the Pauli blocking of electrons is also important as 
T| = 20 K is larger than T = 5 K. For the case of 
T = 200 K, the effect of the Pauli blocking of electrons 
is negligible as T ^> T§. The Pauli blocking of holes be- 
comes visible only for x > 10 -4 , where the hole Fermi 
temperature becomes larger than T — 200 K. 

From Eq. (J5J), one can see that the spin relaxation due 
to the s-d exchange scattering is independent of tempera- 
ture. However, the spin relaxation due to the BAP mech- 
anism increases with increasing temperature because the 
Pauli blocking of electrons and holes decreases with in- 
creasing temperature. Moreover, the matrix element of 
the BAP mechanism increases with the center-of-mass 
momentum of the interacting electron-hole pair, of which 
the ensemble average hence increases with increasing 
temperature>21 The spin relaxation due to the EY mech- 
anism increases with increasing temperature too, as the 



spin-flip matrices [AL k and AL k ] in Eq. ([7]) increase 
with increasing k. Consequently, the BAP and EY mech- 



7 



N Mn ( 10 11 cm" 2 ) 



0.02 0.2 2 20 200 




x 



FIG. 4: (Color online) SRT due to the BAP mechanism 
tbap as function of the Mn concentration at T = 5 (•) and 
200 K (□). Red solid curves: tbap from the full calculation; 
Green dotted curves: tbap from the calculation without the 
Pauli blocking of electrons; Blue dashed curves: tbap from 
the calculation without the Pauli blocking of holes. 



anisms eventually become more efficient than the s-d ex- 
change scattering mechanism at high temperature. 

The appearance of the peak in the t-x curve has been 
observed in a recent experiment at 5 Kj22 However, the 
SRTs we obtain are much larger than the experimen- 
tal value under the same conditions. The deviation may 
come from pretermission of localized holes. At such 
low temperature (5 K), the localization of holes is not 
negligible. 26 ! 62 The localized holes act as exchange in- 
teraction centers located randomly in the sample, which 
thus lead to spin relaxation similar to the s-d exchange 
scattering mechanism. As there is no Pauli blocking 
of the localized holes, the spin relaxation can be very 
efficient^ It should be mentioned that, however, recent 
studies have also shown that there is some compensation 
of the s-d exchange interaction and the electron-hole ex- 
change interaction as holes are always localized on the 
Mn acceptors. 63 ' 64 This leads to a longer spin relaxation 
time^ 3 - and smaller measured (by magneto-optical tech- 
niques) s-d exchange coupling constant 6 ^. However, for 
the high temperature case, the localization is marginal 
and our consideration is close to the genuine case. The 
predicted t-x dependence should be tested experimen- 
tally. 



2. Temperature dependence of the SRT 

We now discuss the temperature dependence of the 
SRT. In Fig. El we plot the SRT as function of Mn con- 
centration x for different temperatures. It is seen that 
for each case the t-x curve shows a peak. It is further 
noted that the temperature dependences of the SRT are 
different for small (e.g., x = 10 -7 ) and large x (e.g., 



x = 10 3 ). To make it more pronounced, we further plot 
the temperature dependence of the SRT for x = 10~ 7 , 
3 x 1CT 5 and 1(T 3 in Fig. [BJ For i = 1(T 7 , the SRT 
first increases then decreases with increasing tempera- 
ture and there is a peak around 20 K. It is understood 
that for such a small x, the electrons and holes are mainly 
from the photo-excitation. For such system, the electron- 
electron and electron-hole Coulomb scatterings are most 
important. It is shown in Ref. [10 that the nonmono- 
tonic temperature dependence of the electron-electron 
Coulomb scattering leads to a peak in the t-T curve. 
In the situation here, the electron-hole Coulomb scatter- 
ing also contributes to the formation of the peak. For 
the case of x = 3 x 10~ 5 , all spin relaxation mechanisms 
are relevant and the most important momentum scatter- 
ing is the electron-impurity scattering. In this case the 
SRT due to the DP mechanism decreases with increasing 
temperature monotonically as the increase of the inho- 
mogeneous broadening dominates* 3 ^, Moreover, the SRTs 
due to the EY and BAP mechanisms also decrease with 
increasing temperature. Consequently, the total SRT de- 
creases with increasing temperature monotonically. For 
the case of large x (x = 10~ 3 ), the spin relaxation is dom- 
inated by the s-d exchange scattering (or the BAP) mech- 
anism at low (or high) temperature. As the s-d exchange 
scattering mechanism is independent of the temperature, 
the temperature dependence is rather weak in the low 
temperature regime. As the temperature increases, the 
EY and BAP mechanisms become more and more im- 
portant, which leads to a fast decrease of the SRT with 
temperature. 

N Mn ( 10 11 cm" 2 ) 
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FIG. 5: (Color online) SRT r as function of Mn concentration 
x at different temperatures. 



3. Photo- excitation density dependence of the SRT 

We now study the photo-excitation density iV cx de- 
pendence of the spin relaxation. In Fig. [7] the SRT is 
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FIG. 6: (Color online) SRT r as function of temperature T 
for different Mn concentrations. 



plotted against the Mn concentration x for three photo- 
excitation densities at low (5 K) and high (200 K) tem- 
peratures. It is noted that the SRT exhibits very dif- 
ferent photo-excitation density dependences at low and 
high temperatures. Moreover, the photo-excitation den- 
sity dependence varies with x. Let us divide the variation 
of x into three regimes: the small x (x < 3 x 10 -6 ) regime 
where the DP mechanism is dominant; the medium x 
(3 x 10~ 6 < x < 10~ 4 ) regime where the DP mechanism 
is comparable with the other mechanisms; the large x 
(x > 10 -4 ) regime where the DP mechanism is irrele- 
vant. 

In small x regime, the DP mechanism dominates the 
spin relaxation. The photo-excitation dependence of the 
DP spin relaxation is different in the degenerate and 
non-degenerate regimes. Similar to the case of n-type 
Ga(Mn)As quantum wells (see Sec. IIIA), in degener- 
ate (low temperature) regime, the density dependence of 
the SRT is dominated by the increase of the inhomoge- 
neous broadening with increasing density, and hence the 
SRT decreases with increasing density. In non-degenerate 
(high temperature) regime, the density dependence of the 
SRT is dominated by the increase of the electron-electron 
and electron-hole Coulomb scatterings with density, and 
hence the SRT increases. 

In the large x regime, the spin relaxation is mainly due 
to the EY, BAP and s-d exchange scattering mechanisms. 
At low temperature, the s-d exchange scattering mecha- 
nism is dominant. As r s d is independent of the electron 
density, the photo-excitation density dependence of the 
SRT (which mainly comes from the EY mechanism) is 
weak. At high temperature, the BAP mechanism domi- 
nates. As holes are mainly from the Mn dopants and the 
electron system is non-degenerate, the SRT also changes 
little with photo-excitation density. 

In the medium x regime, all the four mechanisms con- 
tribute to the spin relaxation. As r s d is independent of 
the electron density, the density dependence comes from 



the other three mechanisms. At low temperature, be- 
sides the DP mechanism, the BAP mechanism is also 
important. However, the BAP spin relaxation changes 
slowly with electron (hole) density as the Pauli block- 
ing is important at low temperature (see Fig. [4J. The 
spin relaxation due to the EY mechanism also increases 

with increasing iV GX as the spin-flip matrices [AlVj. and 

A^, k ] in Eq. (|7|) increase with increasing k. However, 
the EY mechanism is usually less efficient than the DP 
mechanism for x < 10~ 4 . As the DP spin relaxation 
increases with increasing photo-excitation density while 
other relevant mechanisms change slowly or less impor- 
tant than it, the peak moves to larger x with increas- 
ing photo-excitation density as indicated in Fig. EJa). 
Moreover, the total SRT decreases with increasing photo- 
excitation density. At high temperature, as the electron 
system is non-degenerate, the inhomogeneous broaden- 
ing changes slowly with iV ex . However, the screening 
increases with increasing N eK as the carrier density in- 
creases. Hence the momentum scattering (mainly from 
the electron-impurity scattering) decreases with increas- 
ing N ex . Therefore, the SRT due to the DP mechanism 
decreases with increasing N cx . The EY mechanism is 
less important than the DP mechanism in this regime 
[see Fig. EJc)]. Moreover, as the hole system is non- 
degenerate, the spin relaxation due to the BAP mech- 
anism increases with hole density. Therefore, the SRT 
also decreases with increasing photo-excitation density. 



4- Effect of magnetic field on the SRT 

We now study the effect of magnetic field on the SRT. 
The magnetic field is applied parallel to the quantum 
well plane, which is perpendicular to the initial elec- 
tron spin polarization (the Voigt configuration). In tra- 
ditional non-magnetic n-type quantum wells, where the 
electron spin relaxation is dominated by the DP mech- 
anism, the magnetic field in the Voigt configuration has 
dual effects on spin relaxation: (i) elongating the spin 
lifetime by a factor of [1 + (o>lt p ) 2 ] (u>l is the Larmor fre- 
quency, t p is the momentum scattering time) ;— (ii) mix- 
ing the in-plane and out-of-plane spin relaxations 
e.g., ! = + i) when lu l > i(i - Usually, 

effect (i) is weak, but effect (ii) is more important. Differ- 
ing from the case of non-magnetic n-type quantum wells, 
there are several new scenarios in the p-type Ga(Mn)As 
quantum wells: (i) The magnetic field can polarize the 
Mn spins, which alters the spin relaxation due to the s-d 
exchange scattering mechanism, (ii) The nonequilibrium 
Mn spin polarization can be induced during the evolu- 
tion through the s(p)-d exchange interaction with both 
electrons and holes. It precesses around the magnetic 
field and produces the Mn beats. This has been studied 
both experimentally and theoretically in II- VI magneti- 
cally doped quantum wells j 66 ' 67 However, we find that the 
induced nonequilibrium Mn spin polarization is rather 
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FIG. 7: (Color online) SRT r as function of the Mn concen- 
tration for different photo-excitation densities at (a) T — 5 
and (b) 200 K. Red solid curve with •: 7V CX = 0.1 x 10 11 cm" 2 ; 



Green dotted curve with □: N c: 
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dashed curve with A: iV cx = 1 X 10 cm 



small (much smaller than the electron spin polarization) 
and affects the spin dynamics marginally. This is consis- 
tent with the fact that the Mn beats are not observed in 
p-type Ga(Mn)As quantum wells^ (iii) the spin relax- 
ation due to the EY mechanism is also anisotropic be- 
cause the out-of-plane spin relaxation comes from A k 7 



k',k 



while the in-plane relaxation from A k , k [see Eq. 0]. 
At large x when the EY mechanism is important, this 
anisotropy may show up. 

In Fig. we plot the SRT as function of the magnetic 
field with different Mn concentration at T = 5 and 200 K. 
For the case of small x (x = 10~ 6 ), it is seen that the 
SRT increases abruptly when the magnetic field varies 
from to 0.2 T and is almost a constant for B = 0.2 
to 6 T. This abrupt increase of the SRT originates from 
the mixing of the out-of-plane and in-plane electron spin 
relaxations in the presence of magnetic field. For small x, 
the spin relaxation is dominated by the DP mechanism. 
For the DP spin relaxation, the in-plane spin relaxation 
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FIG. 8: (Color online) SRT r vs. the magnetic field with 
different Mn concentrations at (a) T — 5 K and (b) 200 K. 
Note that the values in the figure have been rescaled by a 
factor of 0.25 for the case of x = 10~ 3 at T = 5 K, and 0.05 
(0.4) for the case of x = 10" 4 at T = 5 (200) K. 



is slower than the out-of-plane one, as only part of the 
inhomogeneous spin-orbit field h(k) contributes to the in- 
plane spin relaxation. After the magnetic field is applied, 
the spin relaxation rate becomes i = + ^-). The 

condition for this relation is wr > if- -). In the 

2 v T z Til 1 

situation considered here, it is B > 0.005 T (0.02 T) fol- 
low (high) temperature case. Therefore, the variation of 
the SRT with the magnetic field seems abruptly. 

For the case of large x (x = 10" 3 ), the relevant spin re- 
laxation mechanisms at high temperature are the BAP, 
EY and s-d exchange mechanisms. The BAP and s-d 
exchange mechanisms are isotropic. However, the EY 
mechanism is anisotropic. Our calculation indicates that 

the A^V k term is smaller than the A.^, term in Eq. 0. 
Hence the in-plane EY spin relaxation is faster than the 
out-of-plane one. Therefore the SRT decreases with in- 
creasing magnetic field abruptly by a small amount at 
low magnetic field. After the abrupt decrease, the SRT 
changes little with the magnetic field as the BAP mech- 
anism is almost independent of the magnetic field. At 
low temperature, the dominant spin relaxation mecha- 
nism is the s-d exchange scattering mechanism which is 
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isotropic. Moreover, the contribution of the EY mecha- 
nism is even smaller compared to the high temperature 
case [see Fig- Efa) and (c)]. Therefore, the magnitude of 
the abrupt decrease of the SRT at low magnetic field is 
even smaller than the high temperature case. Let us now 
turn to the magnetic field dependence of the s-d exchange 
scattering mechanism. We choose the eigenstates of a x 
(denoted as |±)) as the basis [hence (s z ) = — J2k^ m Pk\- 
By keeping only the diagonal element of the Mn spin 
density matrix pMm from Eq. ([5]) we obtain 



dt^Pk 

k 



= -JVi 



k 

(S 2 X )] - $><J(e* - e k ){S x )(h, + / k ,_)}. (8) 



As (/k'+ — /k'-) corresponds to the electron spin po- 
larization along the x-axis which is much smaller than 
the Mn spin polarization as both the spin and the g- 
factor of the Mn ions are larger than those of electrons. 
Therefore, the second term in the right hand side of the 
above equation is much smaller than the first one. The 
spin relaxation due to the s-d exchange scattering mech- 
anism increases with increasing magnetic field as (S 12 ) 
does. Consequently, after the abrupt decrease of the SRT 
at low magnetic field, the SRT further decreases with in- 
creasing magnetic field due to the enhancement of the 
s-d exchange scattering. 

We now turn to the medium x case (x — 10 -4 ). At 
low temperature, all the mechanisms are relevant [see 
Fig. |3f a)] . As the BAP and s-d exchange scattering 
mechanisms are isotropic, the anisotropy mainly comes 
from the DP and EY mechanisms. However, as the EY 
mechanism is less efficient than the DP mechanism, the 
anisotropy mainly comes from the DP one. Consequently, 
the SRT first increases abruptly due to the mixing of the 
in-plane and out-of-plane DP spin relaxations, and then 
decreases as the s-d exchange scattering increases with 
increasing magnetic field. At high temperature, the s-d 
exchange scattering mechanism is negligible. Hence af- 
ter the abrupt increase, the SRT varies little with the 
magnetic field. 



IV. CONCLUSION 

In summary, we have performed a systematic investi- 
gation on the spin relaxation in paramagnetic Ga(Mn)As 
quantum wells from a fully microscopic KSBE approach 
with all the relevant scatterings explicitly included. 

For n-type Ga(Mn)As quantum wells, where most Mn 
ions take the interstitial positions ) 23 ! 24 ' 33 we find that the 
spin relaxation is always dominated by the DP mecha- 
nism in the metallic regime. Interestingly, the Mn con- 
centration dependence of the SRT is nonmonotonic and 
exhibits a peak. This is due to the fact that the momen- 
tum scattering and the inhomogeneous broadening have 
different density dependences in the non-degenerate and 



degenerate regimes. A similar effect was found in bulk 
III-V semiconductors very recentlyi 3 - 2 - Our results also are 
consistent with the latest experimental finding that in the 
low Mn concentration regime the SRT increases with Mn 
concentration! 23 ' 24 

For the p-type Ga(Mn)As quantum wells, we study the 
SRT for various Mn concentrations, temperatures, photo- 
excitation densities and magnetic fields. It is found that 
the SRT first increases then decreases with increasing 
Mn concentration. The underlying physics is as follows: 
In the regime of small Mn concentration x, the spin re- 
laxation is dominated by the DP mechanism which de- 
creases with increasing impurity (Mn) density (hence x) 
due to motional narrowing; In the large x regime, as the 
Mn and hole densities are very large, the spin-flip scat- 
terings such as the EY mechanism associated with the 
electron-impurity scattering, the s-d exchange scattering 
and the electron-hole exchange scattering become more 
important than the DP spin relaxation. The SRT hence 
decreases with increasing Mn concentration x and the 
peak is formed. It is found that the most important spin 
relaxation mechanism at large x is the s-d exchange scat- 
tering (or the BAP) mechanism at low (or high) temper- 
atures. The EY mechanism also contributes to the spin 
relaxation at intermediate temperature. 

We also study the temperature dependence of the spin 
relaxation. The behavior also depends on the Mn concen- 
tration x as the relevant spin relaxation mechanisms are 
different for different x. In the small x regime, the SRT 
first increases then decreases with increasing tempera- 
ture which resembles what was found in n-type quantum 
wells with low impurity density^ In the large x regime, 
at low temperature the s-d exchange scattering mecha- 
nism is dominant, which, however, is independent of the 
temperature. The temperature dependence is hence very 
weak. As the temperature increases, the EY and BAP 
mechanisms become more and more important, which 
lead to a fast decrease of the SRT with temperature. In 
the medium x regime, the DP mechanism is also impor- 
tant. As the momentum scattering is dominated by the 
electron-impurity scattering which changes slowly with 
temperature, the increase of the inhomogeneous broad- 
ening leads to the decrease of the SRT. The SRT due 
to the BAP and EY mechanisms also decreases with in- 
creasing temperature. Consequently, the SRT also de- 
creases monotonically with increasing temperature in the 
medium x regime. 

We then address the photo-excitation density depen- 
dence of the SRT. The behavior is different for different 
temperature and x. At low temperature, as the elec- 
tron system is in the degenerate regime, the DP mech- 
anism is largely enhanced as the inhomogenous broad- 
ening increases. However, the s-d exchange scattering 
mechanism is independent of the photo-excitation, and 
the BAP mechanism changes slowly with the photo- 
excitation density as the hole Pauli blocking is very 
strong. The EY mechanism is usually less efficient than 
the DP mechanism. Consequently, the peak in the r- 
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x curve moves to larger x value. The SRTs in the 
small and medium x regimes decrease with increasing 
photo-excitation density as the DP spin relaxation in- 
creases. However, the SRT at large x regime changes 
little as the spin relaxation is dominated by the s-d ex- 
change scattering mechanism. The behavior is quite dif- 
ferent at high temperature where the electron system is 
in the non-degenerate regime. In the small x regime, 
where the momentum scattering is dominated by the 
carrier-carrier Coulomb scattering as the impurity den- 
sity is low. The SRT increases with increasing photo- 
excitation density as the carrier-carrier Coulomb scat- 
tering increases with increasing carrier density. In the 
medium x regime, where the electron-impurity scatter- 
ing is dominant. As the screening increases with increas- 
ing carrier density, the electron-impurity scattering de- 
creases with increasing photo-excitation density. Hence 
the SRT due to the DP mechanism decreases with in- 
creasing photo-excitation density. The EY mechanism is 
less important than the DP mechanism in this regime. 
Moreover, as the BAP mechanism becomes important, 
the SRT decreases with increasing photo-excitation den- 
sity as the hole density increases. In the large x regime, 
the holes mainly come from the Mn dopants and the spin 
relaxation is dominated by the BAP mechanism, the SRT 
hence changes little with photo-excitation. 

We also discuss the magnetic field dependence of the 
SRT. We find that the main effect at low magnetic field 
is the mixture of the in-plane and out-of-plane spin re- 
laxations. The spin relaxation due to the BAP and s- 
d exchange scattering mechanisms is isotropic, whereas 
that due to the DP and EY mechanism is anisotropic. 
For the DP mechanism the in-plane spin relaxation is 



slower than the out-of-plane one, whereas for the EY 
mechanism, the in-plane one is faster than the out-of- 
plane one. Therefore, in small and medium x regimes 
where the DP mechanism is more important than the 
EY mechanism, the magnetic field induces an abrupt in- 
crease of the SRT due to the mixing of the in-plane and 
out-of-plane spin relaxations.— In large x regime, the EY 
mechanism is more important than the DP mechanism, 
and the SRT hence decreases abruptly with increasing 
magnetic field. Another important effect of the magnetic 
field is that it induces an equilibrium Mn spin polar- 
ization and thus enhances the s-d exchange scattering 
mechanism. Consequently, for large x at low tempera- 
ture, where the s-d exchange scattering dominates the 
spin relaxation, the SRT decreases with increasing mag- 
netic field. In other regimes, the s-d exchange scattering 
mechanism is unimportant and the magnetic field depen- 
dence of the SRT after the abrupt jump is also weak. We 
find that the nonequilibrium spin polarization transfered 
from the electron system to the Mn system due to the s- 
d exchange interaction is much smaller than the electron 
spin polarization, which is consistent with the fact that 
the Mn beats are not observed in experiments^ 
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